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Contact inhibition is the process by which cells switch from a motile growing state to a passive 
and stabilized state upon touching their neighbors. When two cells touch, an adhesion link is 
created between them by means of transmembrane E-cadherin proteins. Simultaneously, their actin 
filaments stop polymerizing in the direction perpendicular to the membrane and reorganize to create 
an apical belt that colocalizes with the adhesion links. Here, we propose a detailed quantitative 
model of the role of the cytoplasmic /3-catenin and a-catenin proteins in this process, treated as a 
reaction-diffusion system. Upon cell-cell contact, the concentration in a-catenin dimers increases, 
inhibiting actin branching and thereby reducing cellular motility and expansion pressure. This model 
provides a mechanism for contact inhibition that could explain previously unrelated experimental 
findings on the role played by E-cadherin, /3-catenin and a-catenin in the cellular phenotype and in 
tumorigenesis. In particular, we address the effect of a knockout of the adenomatous polyposis coli 
tumor suppressor gene. Potential direct tests of our model are discussed. 


Introduction 

Before the establishment of cell-cell contacts, epithelial 
cells are in a motile and growing state. The polymerizing 
actin filaments create forces on the membrane that are 
responsible for the formation of lamellipodia and filopo- 
dia [1, 2]. Moreover, the actin filaments undergo contin¬ 
uous branching and growth resulting in dynamic exten¬ 
sions of the membrane [3] . When cells are scarce and do 
not contact each other, E-cadherins are found both on 
the plasma membrane and in membrane vesicles within 
the cytoplasm, but their role is minimal: when located 
on the membrane, they quickly get endocytosed into cy¬ 
toplasmic vesicles [4]. After they have grown enough to 
cover the substrate in a confluent layer, epithelial cells 
become polarized perpendicular to the substrate. At this 
point, they no longer produce lamellipodia and filopodia, 
but instead reorganize their actin into a belt located near 
their apical side (see Fig. 1 A) [5]. Simultaneously, the E- 
cadherins located in the plasma membrane link their ex¬ 
tracellular domains with the cadherins of the neighboring 
cells and colocalize with the actin belt, forming what is 
known as the adhesion zone. The linkage of E-cadherins 
stabilizes their localization on the plasma membrane, ef¬ 
fectively depleting them from the cytoplasm [ 1, 6] . 

The reorientation of the actin filaments upon cell-cell 
contact indicates a reduced activity of branching pro¬ 
teins such as Actin-related proteins 2 and 3 (Arp2/3) 
and an increased activity of bundling proteins such as 
a-catenin dimers (see Fig. 1 B). When oriented parallel 
to it, the growing actin filaments no longer exert a force 
on the plasma membrane. Therefore, the cell downregu- 
lates both its motility and expansion pressure in response 
to reaching confluence, a process referred to as contact 
inhibition. 


In 2005, Drees et al. [7] challenged the textbook view 
that a-catenin mechanically links the adhesion complex 
to the underlying actin cytoskeleton. They showed that 
a-catenin exists either as a monomer or as a dimer, and 
that the domain on an a-catenin monomer that binds 
to /3-catenin and the one that binds to another a-catenin 
monomer overlap. Therefore, the formations of a-catenin 
dimers and a-catenin-/3-catenin complexes are mutually 
exclusive [ ]. Dimeric a-catenin can bundle actin fila¬ 
ments and competes for actin binding sites with Arp2/3. 
According to these findings, a high concentration of 
a-catenin dimers therefore suppresses actin branching, 
growth and expansion pressure (see reviews [9, 10] and 
Fig. 1 B). 

A loss of contact-inhibition via epithelial-mesenchymal 
transition is an essential step for tumorigenesis [11]. It 
has recently been proposed that an excess expansion 
pressure could be a characteristic trait of neoplastic tis¬ 
sues [2]. A breakdown of the regulation mechanism 
discussed above might therefore lead to tumorigenesis. 
It is indeed well-known that the E-cadherin-/3-catenin-a- 
catenin adhesion complex plays an important role in car¬ 
cinomas [13, 14]. A reduced expression of E-cadherins— 
for example due to DNA hypermethylation—is associ¬ 
ated with a loss of cellular polarity and the acquisi¬ 
tion of invasive characteristics [15]. However, it has 
been shown that overexpression and reduced degrada¬ 
tion of /3-catenins also leads to cellular transformations 
that result in the cell’s ability to grow in soft agarose 
gels and to overproliferate at high cell densities [16]. 
Along the same lines, in cells that have undergone the 
epithelial-mesenchymal transition, E-cadherin expression 
is downregulated while the production of /3-catenin is in¬ 
creased [T ]. It has also been shown that the growth- 
inhibiting activity of E-cadherin is counteracted by an 
increased /3-catenin activity [11 ]. Finally, the concern 


2 



FIG. 1: (A) Schematic illustration of the establishment of the 
epithelial cell-cell adhesion zone. After cells have spread via 
protrusions along the substrate and become confluent, they 
start growing upwards and colocalize their actin belt while 
forming the adhesion zone. (B) Different organizations of cor¬ 
tical actin, a-catenin and /3-catenin-related complexes in ep¬ 
ithelial cells during their spreading (left) and after the mature 
epithelial sheet has been formed (right). Before cell-cell con¬ 
tact, /3-catenin-E-cadherin complexes are present in the cy¬ 
toplasm and therefore recruit a-catenin proteins before they 
can form dimers, which lets Arp2/3 complexes branch the 
actin network. In the presence of a neighboring cell however, 
/3-catenin-E-cadherin complexes are mostly found at the cell 
membrane, which favors the formation of a-catenin dimers in 
the cytoplasm. These dimers further bind strongly to actin, 
effectively excluding Arp2/3 complexes from the actin net¬ 
work and favoring parallel bundling. 


tration of /3-catenins is regulated by the adenomatous 
polyposis coli (APC) protein, a tumor suppressor pro¬ 
tein that is known to label /3-catenin for degradation [19]. 
On the other hand, /3-catenin-null cells show an unaf¬ 
fected or even decreased rate of expansion and prolifer¬ 
ation [20, 21]. As the important role of E-cadherin and 
/3-catenin in the progression of cancer has been well stud¬ 
ied, several papers report that the loss of a-catenin is an 
important prognostic factor for cancer, as reviewed in 
Benjamin and Nelson [22]. For example, the ablation of 
a-catenin in the skin causes cellular hyperproliferation, 
occurrence of mitoses away from the basal layer and de¬ 
fects in epithelial polarity [23]. These phenotypes are 
remarkably similar to those obtained with a modified ex¬ 
pression level of E-cadherin or /3-catenin proteins. 

Although strong indications exist that the influence of 


a-catenin on actin polymerization plays a role in con¬ 
tact inhibition, the functional details of this mechanism 
remain unclear. Important progress has been made in 
this direction by Drees and coworkers [7, 8], who pro¬ 
posed a picture in which cell-cell contact leads to an ac¬ 
cumulation of E-cadherin-/3-catenin-a-catenin complexes 
at the adhesion sites. They propose that the release of 
a-catenin monomers from these complexes into the cyto¬ 
plasm provide an increase in a-catenin dimer concentra¬ 
tion, favoring actin bundling and downregulating actin 
assembly and branching. In the present work, we propose 
a model for the E-cadherin-/3-catenin-a-catenin function 
that is based on a reaction-diffusion system. We show 
that the interplay between these three proteins results 
in a pathway for contact inhibition that downregulates 
actin polymerization in response to cell-cell contact. 

Our mechanism relies on the fact that the binding of 
/3-catenin to a-catenin limits a-catenin dimerization in 
the cytosol. When /3-catenin-E-cadherin complexes are 
recruited to the cell membrane due to cell-cell contact, 
the cytosolic concentration of /3-catenin drops and a- 
catenin dimerization can take place. According to the 
work by the Nelson and Weis group [7], this in turn pre¬ 
vents Arp2/3-based actin branching and cause the cell to 
enter a quiescent state. Using the framework of our phys¬ 
ical model, we investigate the effect of disruptions of this 
pathway and obtain results that are consistent with ex¬ 
perimentally observed cellular transformations that lead 
to tumorigenesis. 


Results 

Description 

The mechanism of the pathway we propose is schemat¬ 
ically illustrated in Fig. 2. Its main feature is that 
a-catenin-/3-catenin binding competes with a-catenin 
dimerization: At high cytosolic concentrations of /3- 
catenins, the majority of a-catenins enters a-catenin-/3- 
catenin complexes. At low cytosolic concentrations of 
/3-catenins however, a-catenins form dimers almost ex¬ 
clusively [7]. Therefore, the organization and activity of 
the actin cortex of the cell depends on the presence of a 
neighboring cell by the following mechanism: It is known 
that /3-catenin quickly binds to E-cadherin after produc¬ 
tion at the Golgi apparatus of the cell [24] . When the cell 
is in its growth phase, E-cadherin-/3-catenin complexes 
are mostly found in vesicles in the cytoplasm [4], effec¬ 
tively creating a large concentration of /3-catenin com¬ 
plexes in the cytosol. These complexes further recruit 
most of the a-catenin monomers that are present in the 
cytoplasm, letting actin binding sites free for Arp2/3 
complexes to bind. The structure of the actin cortex 
is therefore branched, and its activity is high. In con¬ 
trast, when contact with a neighboring cell is established, 
the E-cadherins bind to the neighboring cell and accu¬ 
mulate at the membrane, effectively lowering their con- 
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FIG. 2: Schematic illustration of the proposed cadherin- 
catenin pathway for contact inhibition as well as its possible 
disruptions. Arrows and T-bars between the different gene, 
proteins or cell states of the diagram indicate induction and 
repression, respectively. The signs of different colors illustrate 
how various events can lead to a breakdown of this path¬ 
way: reduced expression of E-cadherins (red, upper signs), 
mutation of the APC tumor suppressor gene ( orange, middle 
signs), and reduced expression or mutations of a-catenin ( yel¬ 
low, lower signs). Minus signs indicate either decreased con¬ 
centrations or complete impairment of the associated proteins, 
and plus signs indicate increased concentrations as compared 
with the healthy cell state. 


centration in the cytosol [25]. Since a large fraction of 
E-cadherins is bound to /3-catenins, the establishment 
of cell-cell contacts also induces a redistribution of /3- 
catenins to the plasma membrane. Indeed, the potent 
ability of E-cadherin to recruit /3-catenin to the cell mem¬ 
brane has been observed in vivo [26]. This /3-catenin re¬ 
distribution to the plasma membrane in turn favors the 
formation of a-catenin dimers in the cytosol, which fur¬ 
ther favors actin bundling rather than actin branching 
and polymerization. Note that other protein complexes 
could play a similar role as E-cadherins and transport 
/3-catenin proteins to the cell membrane upon cell-cell 
contact, as it has been proposed recently [27]. In any 
case, the cell switches from an active state with high 
actin branching and polymerization activity when there 
is no cell-cell contact, to a passive state characterized by 
a reduced actin activity after cell-cell contact has been 
established. 


Model Equations 

To model the cadherin-catenin system described above 
in a quantitative manner, we write a system of reaction- 
diffusion equations for the cytosolic concentrations of the 
different proteins involved. (Note that active transport 
of these proteins may be involved but that we do not 
expect the mechanism presented in this paper to depend 
crucially on this aspect.) In this model, we treat all pro¬ 
tein bindings as irreversible because binding affinities are 
high (typically with energies of many k^T [28]). Assum¬ 
ing that the protein production rates and the configu¬ 
ration of neighboring cells are constant in time, we can 
focus on the steady-state dynamics of the system, which 
for the cytosol of the cell can be written as 

D a \7 2 C a -k a pC a Cp-2k aa C 2 a -r a C a = 0 (la) 
DpV 2 Cp-k a pC a Cp-rpCp = 0 (lb) 

C (xol + kaoiC^ (t’o'O: H - ^aa) ^aa — 0 (lc) 

Cq,( 3 k a pC a Cp 't'afiCoLp = 0. (Id) 

These equations respectively describe the diffusion dy¬ 
namics of a-catenin, /3-catenin, a-catenin dimers and 
a-catenin-/3-catenin complexes in the cytoplasm. Here 
C a , Cp, C aa and C a p are the cytoplasmic protein con¬ 
centrations of a-catenin monomers, /3-catenin monomers 
(bound to cytosolic E-cadherins), a-catenin dimers, and 
a-catenin-/3-catenin complexes (bound to cytosolic E- 
cadherins), respectively; D a , Dp, D aa and D af 3 are the 
associated diffusion constants and r a , rp, r aoi and r a p 
the associated degradation rates; k aa and k a p are re¬ 
spectively the rates of a-catenin dimerization and a- 
catenin-to-/3-catenin binding; and f aa is the reaction rate 
of a-catenin dimers with actin. Note that since most 
/3-catenins bind to E-cadherins immediately after pro¬ 
duction [2 ], we do not explicitly model the reaction- 
diffusion dynamics of E-cadherins but instead account for 
its important effect on the redistribution of /3-catenin-E- 
cadherin complexes in the boundary conditions at the 
plasma membrane, as we shall see below. Modeling the 
diffusion and reactions of E-cadherins and /3-catenins sep¬ 
arately adds another layer of complexity to our model, 
but would not qualitatively change our main results. 
Therefore, when we refer to /3-catenin in our model, we 
implicitly mean the E-cadherin-/3-catenin complex. Note 
also that the effect of the Wnt signaling pathway on /3- 
catenin is taken into account effectively in the bulk degra¬ 
dation rate of this protein. 

Production of these proteins in the vicinity of the cell 
nucleus as well as their interactions with the plasma 
membrane need to be accounted for using appropriate 
boundary conditions. The production of a-catenin and 
/3-catenin in the Golgi apparatus of the cell is taken 
into account by fixed influxes of proteins into the cyto¬ 
plasm, denoted by and jp, respectively. On the mem¬ 
brane, the concentrations of protein complexes are (7™’ d , 
C™’ a , C™g d , C™g a —all bound to E-cadherin proteins— 
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which can be either detached or attached to an adja¬ 
cent cell via E-cadherin-E-cadherin homophilic binding 
as indicated by the superscripts d and a. Cytoplasmic 
concentrations at the membrane, denoted Cjj, and 
Cl p , correspond to the respective concentrations intro¬ 
duced in Eq. 1 at this particular location. /3-catenin as 
well as a-catenin-/3-catenin complexes—both bound to 
E-cadherin proteins—can go to the plasma membrane of 
the cell where they are then in the detached state. We 
denote by fc™ and kp S ( k °(g and fc°®) the rates at which 
the protein complex /3-catenin-E-cadherin (a-catenin-/3- 
catenin-E-cadherin) goes to the plasma membrane of the 
cell, and note that only complexes in the detached state 
can move from the membrane to the cytoplasm. The two 
fluxes, j/ 3 and j a p , of /3-catenin and a-catenin-/3-catenin 
complexes from the cytoplasm to the plasma membrane 
of the cell then read 

3P = k°fC*-kfCf’\ (2a) 

lap = kZfrC'p - kfpCZf. (2b) 

Monomeric a-catenins can only go to the plasma mem¬ 
brane by forming a-catenin-/3-catenin complexes via a 
reaction with /3-catenin complexes that are already lo¬ 
cated on the membrane (either in the attached or de¬ 
tached state). The rates of these reactions are denoted 
respectively by fc™^ a and . Once formed, these 
complexes do not release pure a-catenins anymore, a- 
catenin dimeric complexes on the other hand cannot 
go to the membrane since they neither attach directly 
to E-cadherins, nor can they bind /3-catenin-E-cadlierin 
complexes because for that they need to be in their 
monomeric form. Their flux therefore vanishes. We fi¬ 
nally get 

Ja — C 'p "V K oc )3 °/3 u o) 

j aa = 0. (3b) 

To solve our system of equations, we must combine 
these boundary conditions with the cytosolic protein dif¬ 
fusion equations (Eqs. la-ld), which we do thanks to 
the definition of diffusive fluxes {j a = —Da^Ca)- To do 
so, we need to eliminate the membrane protein concen¬ 
trations from our system of equations. This is done by 
writing the balance of protein complexes located on the 
plasma membrane of the cell: In addition to the reac¬ 
tion rates introduced above, we introduce the rate 
for a given complex A linked to an E-cadherin molecule 
to attach to another complex B of the adjacent cell via 
cross-membrane E-cadherin homophilic binding. Also, 
all complexes are degraded with their specific rates r™’ d 
and r™’ a on the membrane. For simplicity, we assume 
a completely symmetric, identical configuration of the 
neighboring cell, and thus identical protein concentra¬ 
tions on the membrane of the adjacent cell. Taking all 
of this into account, the protein concentrations on the 
cell membrane are determined by the following steady- 
state equations, a schematic representation of which is 


presented in Fig. 3: 
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Here, the two first equations describe the balance of 



FIG. 3: Schematic illustration of the reactions occurring at 
the cell membrane and leading to Eq. 4. E-cadherin vesicles, 
either bound to /3-catenins alone or to a-catenin-/3-catenin 
complexes, can merge with the membrane or be endocy- 
tosed. The two associated light gray arrows correspond to 
the four on and off rates in Eq. 4. In the presence of cell¬ 
cell contact, the E-cadherins on the membrane can bind to 
E-cadherins on the membrane of the adjacent cell, which is 
represented by the black arrows for /3-catenin-associated com¬ 
plexes, and by the dark gray arrows for a-catenin-/3-catenin- 
associated complexes. These correspond to all the reactions 
which have rates with an EE superscript in Eq. 4. In addi¬ 
tion, a-catenin monomers can bind to E-cadherin-/3-catenin 
complexes present on the membrane, whether they are bound 
or not to E-cadherins from the adjacent cell. This is repre¬ 
sented by the vertical arrows and corresponds to the terms n 
Eq. 4 whose rates have an m,d or m,a superscript. Finally, all 
protein complexes located on the membrane can be degraded, 
which is taken into account by the rates labeled with the let¬ 
ter r in Eq. 4 (not represented). We assume a symmetric 
configuration of the adjacent cell. 

/3-catenin and a-catenin-/3-catenin complexes, respec¬ 
tively, on the plasma membrane of the cell that are de¬ 
tached from the neighboring cell, and the two follow- 
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ing equations do the same for the attached protein com¬ 
plexes. For example, in the first equation, /3-catenin- 
E-cadherin complexes in the detached state can- in the 
order of the terms present in the equation- be replen¬ 
ished via attachment of /3-catenin-E-cadherin complexes 
from the cytoplasm, disappear via endocytosis of /3- 
catenin-E-cadherin complexes, form a-catenin-/3-catenin- 
E-cadherin complexes, attach with other /3-catenin-E- 
cadherin complexes from the neighboring cell via cross¬ 
membrane E-cadherin-E-cadherin homophilic binding, or 
with a-catenin-/3-catenin-E-cadherin complexes from the 
neighboring cell, and, finally, disappear via degradation. 
The terms in the second equation are of similar origin. 
The third and fourth equations for the attached states 
resemble the previous ones, except that there is no ex¬ 
change of proteins directly with the cytoplasm in this 
case. 


Steady-State Concentration Profiles 


To solve the system of equations given by Eqs. 1-4, we 
now separate the two cases of presence and absence of 
contact with a neighboring cell, for which we can sepa¬ 
rately eliminate the membrane concentrations from the 
boundary conditions given by Eqs. 2 and 3 thanks to 
Eq. 4. In the absence of cell-cell contact, the different 
rates fc EE s vanish and all the proteins on the membrane 
are in the detached state. In this case, Eqs. 4c and 4d 
become irrelevant, and Eqs. 4a and 4b become 
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The first equation allows to solve for C™’ d and then ex¬ 
press j a and jp as a function of the cytosolic concentra¬ 
tions of a-catenin and /3-catenin complexes only. Thus, 
we get a closed set of equations for these two quantities, 
in which we find the cytosolic equations (Eqs. la and lb), 
as well as the expressions for the fluxes at the boundaries 
of the system, namely the imposed fluxes and jp at the 
Golgi apparatus of the cell and the two following fluxes 
at the cell membrane, 
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Jet — 
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(6a) 
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This system can be solved independently and then used 
to solve for the concentrations of a-catenin-/3-catenin and 


a-catenin dimeric complexes in a second step, using the 
remaining equations. 

In the presence of cell-cell contact, we assume for sim¬ 
plicity k EE —► oo, meaning that all protein complexes on 
the membrane instantaneously bind to the neighboring 
cell. Therefore, the concentrations of unbound proteins 
on the membrane C'S 1 ’ d and C™p d vanish. Similarly to the 
previous case, the dynamics for a-catenin and /3-catenin 
complexes decouples from the rest of the system, and the 
fluxes can be obtained from Eqs. 2 and 3 after we have 
solved for C“’ a using Eqs. 4a and 4c: 
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The system of equations derived above can now be 
solved independently in the two configurations of the 
cell numerically, namely in the presence or absence of 
contact with a neighboring cell. It consists of Eqs. la 
and lb, together with the boundary conditions (Eq. 6 
or 7) at the plasma membrane and a constant protein 
influx given by and jg at the Golgi apparatus of the 
cell. We solve this system for the case of a spherical 
cell of radius R, whose Golgi apparatus is modeled as a 
sphere of radius ro- In Fig. 4, we illustrate the differ¬ 
ence in the concentrations of a-catenin, /3-catenin and 
a-catenin dimers with and without cell-cell contact. We 
see in Fig. 4 C that the overall concentration of a-catenin 
dimers presents a significant increase in the case of cell¬ 
cell contact as compared to the case without contact, 
which provides an efficient switch between the two phe¬ 
notypic states of the cell. Note that in both cases, there 
is a drop in the concentration of a-dimers away from the 
nucleus. If the diffusion constant is small enough (i.e. 
if D aoi /[R 2 (r aa + r aa )] <C 1), this concentration drop is 
significant and could be relevant for the spatial organiza¬ 
tion of polymerized actin within the cell. We comment 
further on this aspect in the discussion section. 


Scaling Analysis 

Let us now perform a scaling analysis of the total num¬ 
ber of a-catenin dimers in the system as given by our 
model, comparing the two cases with and without cell¬ 
cell contact. In asymptotic limits in which the involved 
length scales separate, it is possible to solve our system 
of equations analytically. We thereby obtain a better 
physical understanding of the contact inhibition mech¬ 
anism proposed in this paper. We also derive a sim¬ 
ple expression for the change in the total amount of a- 
catenin dimers N aa in the cell between the contact and 
no-contact states, which dictates the amplitude of the 
switch. This final expression is given by Eq. 19, and 
one may want to skip to this equation and its associated 
comments directly. Later, this derivation also helps us 
to see in which biological conditions our mechanism can 
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FIG. 4: Cytosolic concentration profiles of /3- catenin (A), 
a-catenin (B) and a-catenin dimers (C) resulting from the 
reaction-diffusion system described by Eqs. la-ld, with the 
boundary conditions without cell-cell contact (Eq. 6) or with 
cell-cell contact (Eq. 7), as functions of the distance from the 
center of a cell with spherical symmetry. The different param¬ 
eters are as follows. The Golgi apparatus of the cell is located 
at ro = 1 /im and the total radius of the cell is R = 10 fun. 
The /3-catenin influx is jg = 5.0 /im /dVl s -1 . The diffusion 
constants of a- and /3-catenin are equal to 1 /im 2 s -1 , and the 
one for a-catenin dimers is 0.5 /im 2 s _1 [30]. In these plots, 
the reaction rate of a-catenin with /3-catenin is 0.01 /tM -1 s -1 , 
and the one of a-catenin with itself is 0.005 /jM _1 s _1 [31]. 
The protein degradation rates of a- and /3-catenin are equal 
to 10 -3 s _1 and the consumption and degradation rates of a- 
catenin dimers are 0.5 10 -3 s _1 . Finally, the membrane bind¬ 
ing and unbinding rates fcS n and k° B of the /3-catenin protein 
complex are equal to 1 /im s _1 and Is -1 , respectively. 


function, as well as to investigate the various possibilities 
that can lead to its breakdown. This is done in the next 
section. 

We first identify the different characteristic lengths 
over which the concentrations of the different proteins 
vary as each of the reactions is considered separately. 
For a given protein, the shortest of the characteris¬ 
tic lengths of the different reactions it enters deter¬ 
mines its dominant reaction pathway. The character¬ 


istic length for the change in monomeric a-catenin con¬ 
centration due to a-catenin dimer formation (a-catenin- 
to-/3-catenin binding) is given by Z“ a = yjD a /(k aa C a ) 
(Z“o = yjD a /(k a pCp)). In a similar way, the change 
in /3-catenin concentration due to a-catenin-to-/3-catenin 
binding is given by = \J Dp / (k a pC a ). Finally, the 
characteristic length due to monomeric a-catenin degra¬ 
dation (/3-catenin degradation) is l a = yjD a /r a (Ip = 
\/Dp/rp). 

We first look at the case where is the shortest of 
the length scales given above. As we shall see below, 
such a condition is realized as soon as the production 
of /3-catenin at the Golgi apparatus of the cell is large- 
enough, such that reactions with monomeric /3-catenin 
proteins (bound to E-cadherins) are fast. In this case, a- 
catenin-to-/3-catenin binding is dominant over a-catenin 
dimerization in the absence of cell-cell contact, such that 
our mechanism can be efficient. (Other limits are studied 
below.) The concentration of a-catenin at the cell mem¬ 
brane is always very low and the change in the steady 
state concentration of a-catenin with and without cell¬ 
cell contact results from a redistribution of /3-catenin in¬ 
side the cell. 

Within this limit, we can assume a quasi-constant con¬ 
centration of /3 catenin in a region of length l^p around its 
source, which allows us to find analytical expressions for 
the reaction-diffusion system in a one-dimensional geom¬ 
etry with coordinate x, the protein source being at x = r 0 
and the cell membrane at x = R. The solution for the 
a-catenin concentration is given by: 

C a ~ exp (—(x — r 0 )/lap) , (8) 


with C° = l^pfJDa = yjk a pCpD a , and where C° 

and Cp are the concentrations of a-catenin and /3-catenin 
complexes, respectively, at x = Vq- The solution for the 
/3-catenin concentration is given by 


Cp ~ Cp cosh ((x - ro)/Ip) - 


1 ° 

Jp 


'Jo 


sinh ((x - r 0 )/lp).(9) 


y/Dprp 

From the boundary conditions, we can determine the ex¬ 
pression for Cg: 

„-o 




Jp 


V D P r P ’ 


( 10 ) 


where 

1 + 5 tanh (R/lp) 


f = 


tanh (R/lp) + g 


and g = 


~off 


V D P r P 


( 11 ) 


It is now possible to translate our initial assumptions 
on the different characteristic reaction lengths into con¬ 
ditions directly on the concentration C°p of /3-catenin at 
x = r 0 . For the three characteristic lengths h = l a ,lp 
and R , the condition li reads 


C?» 


Da 1 

kap li' 


( 12 ) 
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while the conditions Pp <C //g and Pp <C Z“ a , respec¬ 
tively, read 


and 


C?» 



C?» 


(2 fcggj°) 2/3 
fc «/3 ( A *) 1/3 


(13) 


(14) 


Finally, there is an additional condition stating that the 
concentration of /3-catenin is quasi-constant in a region 
of length around the protein source: 




D a rp 1 

k a pDp f 1 


(15) 


Since oc jp — , all of these conditions are satisfied 

for a sufficiently large influx j% of /3-catenins at the Golgi 
apparatus of the cell. 

The amplitude of the switch is given by the change 
in the total amount of a-catenin dimers N aa in the cell 
between the contact and no-contact states. The concen¬ 
tration in a-catenin dimers is simply given by 


Caa — 


r 2 

^ a. i 


(16) 


if P y p is much smaller than the two other lengths scales 
given by Eq. lc, namely = \JD aa C aa /{k aa C 2 ) and 
l aa = a JD aa /(r aa + r aa ). Integrating the a-catenin 
dimer concentration over the size of the whole cell un¬ 
der this hypothesis, we obtain 


N„ 


k G 0 ! 2 

2 r aa (k a p) z / 2 (D a ) 1 / 2 


(C°p )- 3/2 


(17) 


as a formal expression. The consistency check for this 
expression gives the following condition on (72: 


Cl » 


/4,. k„ 


Daa P o 


(18) 


which comes from the dependence of N aa on C'p, and 

where /»° ff ’ nocon j s the off-rate of /3-catenin from the mem¬ 
brane when their is no cell-ccll contact. Hence, for a pro¬ 
tein degradation rate r™’ d much smaller than the off rate 
kp S , we expect a significant switch in the total amount of 
a-catenin dimers produced and, thus, a functional con¬ 
tact inhibition mechanism. 

We now look at two cases where Pp may not neces¬ 
sarily be the smallest characteristic length in the system. 
First, it is possible that, when there is cell-cell contact, 
1“ 0 becomes the shortest characteristic length instead of 
Pp in the absence of contact. Indeed, when there is cell¬ 
cell contact, /3-catenin proteins could be sufficiently de¬ 
pleted from the cell for most of the a-catenins to form 
homo-dinreric complexes before reacting with /3-catenins. 
In this case, the contact inhibition switch remains intact, 
and the previous ratio still scales as stated in Eq. 19. 

Another limit corresponds to the case where the cell 
radius, R, is the shortest length scale in the system. In 
this case, the system of reaction diffusion equations— 
together with the corresponding boundary conditions— 
can be treated as a system without spatial extension. 
A substantial change in the /3-catenin concentration be¬ 
tween the two states of the cell can be achieved in this 
limit if the degradation rate of /3-catenin in the cytosol 
rp is much smaller than its degradation rate on the mem¬ 
brane r^' d . We present the numerical solutions that cor¬ 
respond to this limit in Fig. 5. In particular, Fig. 5 B 




0.04 0.06 

kf [1/s] 



which again is satisfied for a sufficiently large influx jp of 
/3-catenins at the Golgi apparatus of the cell. The results 
in the presence and absence of cell-cell contact can be 
obtained by switching the rates for detached membrane 
proteins with those for attached ones. In particular, the 
off-rate of /3-catenin from the membrane, kp S , must be set 
to 0 in the case where there is contact with a neighboring 
cell. If all the other rates stay the same, we obtain a 
simple expression for the ratio of the total amounts of 
a-catenin dimers with and without cell-cell contact: 


jvrcon 

aa 


1\J nocon 
1 y aa 



i off,nocon 



3/2 


(19) 


FIG. 5: Numerical solutions for the total number of a- 

catenin dimers in the cell as a function of different param¬ 
eters when the cadherin-catenin system is treated as a zero¬ 
dimensional reaction system. In this limit, the contact in¬ 
hibition mechanism is based on a larger degradation rate 
on the membrane (rT’ d = 10 -2 s -1 ) than in the cytosol 
(rp = 1CF 2 s _1 ). The a-catenin degradation rate r a is also 
assumed to be small (10~ 2 s -1 ). The total production rates of 
a-catenin and /3-catenin are respectively equal to 7.5 10' 5 s^ 1 
and 14 10 3 sW Protein-protein reaction rates are 10 -3 s _1 , 
and on- and off-rates of the /3-catenin protein complex to 
the cell membrane are kp n = kp S = 0.1 s _1 . Note that for 
rp > r™’ d , the switch is reversed. 

















shows that the switch is controlled by the ratio rp/r™' d . 
As an alternative, the reaction rate of membrane- 
bound /3-catenin with a-catenin could be smaller than 
the corresponding reaction rate in the cytosol k a p. This 
would also yield a functioning contact inhibition switch. 


Breakdown of Contact Inhibition 


Let us now investigate the various possibilities that, 
according to our model, can lead to a breakdown of the 
contact inhibition pathway. To investigate what affects 
the production of a-catenin dimers, let us look at the 
total amount of a-catenin dimers in the cell calculated 
from a numerical solution of our whole system of equa¬ 
tions (Eqs. la-ld) together with the boundary condi¬ 
tions described by Eqs. 6 and 7, for the contact-free and 
contact-inhibited states (see Fig. 6). We first show the 



FIG. 6: Integrated a-catenin dimer concentration over the 
cell volume from the numerical solution of the whole model. 
This quantity is plotted as a function of the a-catenin influx 
ja (A), of the /3-catenin degradation rate in the cytoplasm rp 
(B) and of the unbinding and binding rates of the /3-catenin 
protein complex off and on the membrane kp S (C) and kp n 
(D). The same constants are used as in Fig. 4. In B, r™’ d = rp 
is assumed. Contact inhibition occurs when the concentra¬ 
tion of a-catenin dimers C aa is large. The contact-inhibition 
breaks down for low production j% of a-catenin (A) and for an 
increased degradation of /3-catenin (B). It also breaks down 
for an increased membrane off-rate (C), which corresponds to 
mutations of E-cadherins leading to less efficient formation 
of cell-cell E-cadherin bonds. This results in an insufficient 
trapping of /3-catenins on the membrane in the presence of 
cell-cell contact. Finally, contact inhibition breaks down for 
small values of fc™ (D), which could correspond to a less effi¬ 
cient binding of /3-catenin to the membrane due to a decreased 
expression of E-cadherins. 

dependence on (Fig. 6 A), which is proportional to 
the total production of a-catenin in the cell Golgi appa¬ 


ratus. As discussed above, many cancerous cells show 
mutations impairing the function or production of a- 
catenin proteins. From Eq. 17, we see that the total 
amount of a-catenin dimers scales like (j°) 2 , which is 
consistent with these experimental observations. For low 
values of the difference between the contact-free and 
contact-inhibited state disappears, as it has been exper¬ 
imentally observed [23]. Next, we consider the effect of 
a knockout of the APC protein, which is known to la¬ 
bel the /3-catenin in the cytosol for degradation. From 
Eqs. 10, 11, and 17, we see that for a fast /3-catenin 
degradation, the total amount of a-catenin dimers scales 
like (rp) 3 / 4 . Hence, contact inhibition becomes less ef¬ 
fective for a lower cytosolic /3-catenin degradation rate, 
as it has been observed in experiments [19]. As can be 
seen in Fig. 6 B, a low degradation rate of /3-catenin 
which corresponds to a knockout of APC—leads to a to¬ 
tal concentration of a-catenin dimers in the cell that is 
low, even when there is cell-cell contact. Other defects 
frequently encountered for cancerous cells are downregu- 
lation or mutations of E-cadherins [15]. In our picture, a 
malfunction of E-cadherins due to mutations corresponds 
to a less effective binding of E-cadherins to neighboring 
cells, and thereby a less efficient trapping of E-cadherins 
and /3-catenins to the plasma membrane. Fig. 6 C shows 
that an increased off-rate in both states again leads 
to a failure of contact inhibition, since the difference be¬ 
tween the contact and no-contact states disappears for 
large values of kp S in the contact state. Finally, a lower 
expression of E-cadherin results in a less effective bind¬ 
ing of /3-catenin to the plasma membrane, which can be 
modeled by a decreased on-rate k™, as shown in Fig. 6 D. 
Similar results are obtained when R is the shortest length 
scale in the system. In this case, the solution of the sys¬ 
tem of equations is plotted in Fig. 5 as a function of the 
model parameters that simulate a breakdown of the con¬ 
tact inhibition mechanism. The breakdown of the switch 
is similar to the one discussed above. 


Discussion 

In this article, we have proposed a reaction-diffusion 
model of the cadherin-catenin system in which the con¬ 
centration of a-catenin dimers increases in a confluent 
cell as compared to a cell without contact. We propose 
that this switch is due to a competition between mutu¬ 
ally exclusive a-catenin dimerization and a-catenin-to-/3- 
catenin binding in the cytosol of the cell. In the presence 
of cell-cell contact, intercellular E-cadherin bonds pre¬ 
vent endocytosis of E-cadherin complexes. This leads to 
a redistribution of unbound /3-catenins to the cell mem¬ 
brane and thereby a significant increase in the amount of 
a-catenin dimerization. Hence, the cell shifts between an 
active state combining high cellular expansion pressure 
with high cellular motility, to a quiescent state, where 
actin branching is inhibited. 

From our analysis, we expect the contact inhibition 











9 


switch to function efficiently if /3-catenin is sufficiently 
abundant in the cell to effectively compete with a-catenin 
dimerization. Therefore, there are three distinct possi¬ 
bilities that can lead to a functioning contact inhibition 
mechanism. First, the protein reaction rates could be 
sufficiently fast compared with protein diffusion to effec¬ 
tively separate two distinct pools of /3-catenin proteins, 
respectively cytosolic and membrane-bound /3-catenins 
(both linked to E-cadherin proteins). While /3-catenins 
in the cytosol compete with a-catenin dimerization, /3- 
catenins on the membrane cannot react with a-catenins 
anymore, because all a-catenins either bind /3-catenins 
or form dimers before they get a chance to arrive at the 
membrane. In that case, the contact inhibition switch 
comes from a redistribution of /3-catenin-E-cadherin com¬ 
plexes from the cytosol toward the cell membrane as a 
response to contact with a neighboring cell, effectively 
letting a-catenin dimers form in the cytosol before reach¬ 
ing the cell membrane, where a high concentration of /3- 
catenin proteins is found. Second, the contact-inhibition 
switch could arise from depletion of /3-catenins from the 
entire cell in the state with contact as compared to the 
state without contact. This is case for example if the 
degradation rate of /3-catenin is much larger on the mem¬ 
brane than it is in the cytosol, which is possible since 
/3-catenin degradation takes place via two distinct path¬ 
ways in the cytosol and on the membrane of the cell. In 
this case, the mechanism works even for very slow re¬ 
action rates and fast protein diffusion, and the spatial 
structure of the cell can then be ignored or be treated 
as a zero-dimensional system. This case is illustrated in 
Fig. 5. Third, the switch could arise from a reaction of 
membrane-bound /3-catenin with a-catenin much slower 
than the corresponding reaction rate in the cytosol. De¬ 
pending on the state of the cell, /3-catenin proteins are 
indeed located primarily either on the plasma membrane 
or in the cytosol, allowing a-catenin proteins to dimerize 
or not. This case also does not rely on slow diffusion and 
can be treated as a zero-dimensional system. 

In addition to providing a mechanism for contact in¬ 
hibition, the model qualitatively reproduces the effect of 
several mutations that are known to cause the breakdown 
of this mechanism and result in tumor-like phenotypes. 
As we have seen from Figs. 5 and 6 , the model agrees 
with experimental observations in which the expression 
level or the degradation rate of either E-cadherin, /3- 
catenin or a-catenin are modified. These findings ex¬ 
plain why a broad range of mutations leads to similar 
cancerous phenotypes. In particular, the effects of an 
increased /3-catenin concentration on contact inhibition 
are explained by this model without implicating the Wnt- 
signaling pathway [18]. 


While this work shows that the cadherin-catenin 
reaction-diffusion system could play the role of a contact 
inhibition switch, it is impossible to determine whether it 
is the most relevant effect without further experimental 
studies. Experiments that would directly test this path¬ 
way are possible. We have already discussed how the 
cadherin-catenin mechanism reproduces the observed ef¬ 
fects of a change in the production and degradation rates 
of different proteins. However, it might also be possible to 
only inhibit the interaction of any given pair of these pro¬ 
teins by phosphorylation of specific residues, and thereby 
directly investigate every step in the proposed mechanism 
without interfering with other pathways like the Wnt- 
signaling pathway [18, 29]. For example, one experiment 
of particular interest would be to knockout /3-catenin-to- 
a-catenin binding by phosphorylation without changing 
the level of expression of these proteins. This would dis¬ 
tinguish our model from the picture proposed by Nelson 
et al. [7, 8]: indeed, within the reaction-diffusion model 
presented here, such a treatment would result in an in¬ 
creased concentration of a-catenin dimers, the inhibition 
of actin branching and thus the contact-inhibited cell 
state even in the absence of a neighboring cell. In con¬ 
trast, in the picture proposed by Nelson et al., failure of 
/3-catenin to bind to a-catenin would lead to a disruption 
of the localization of a-catenin to the adhesion sites, and 
thereby to actin branching and polymerization even at 
confluence. 

If the model presented here were to be confirmed exper¬ 
imentally, it could potentially ground the idea that differ¬ 
ent homeostatic growth pressures between neoplastic and 
healthy tissues are responsible for tumor growth [12]. In 
particular, one could then test whether different disrup¬ 
tions of the cadherin-catenin pathway that are known to 
lead to tumorigenesis would affect the homeostatic pres¬ 
sure growths of the tissues under study. Such an obser¬ 
vation could potentially give a direct explanation of the 
observed link between the cadherin-catenin system and 
neoplastic phenotypes. 
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